Two exercises about neutrino departure times at CERN 
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Two simple exercises are solved, which educators can use to awake interest of their students in 
subtleties of the CERN Neutrino beam to Grand Sasso (CNGS) experiment. The first one is about 
the statistical error of the average departure time of neutrinos from CERN. The second one about 
a hypothetical bias in the departure times. 

PACS numbers: 01.40.-d, 29.25.-t 



In a highly publicized measurement [l| neutrinos from 
CERN arrived at Gran Sasso about 



St = 60 ns (nano seconds) 



(1) 



too early. Here we perform two pedagogical exercises 
using numbers given in the paper. 

The neutrinos are produced in extractions that last 
about 10,500 ns each, smashing protons into a graphite 
neutrino production target. About 6.2 x 10 15 protons par- 
ticipate in one extraction and any proton could have pro- 
duced the finally observed neutrino. Altogether 16,111 
neutrino events entered the analysis. 

Exercise 1: Assume the proton distribution is uni- 
form in time, what is the initial uncertainty in the neu- 
trino creation time? 

This is easy, let us center the time interval about zero, 
i.e., from [-t max ,*max] with t max = 5,250ns. The vari- 
ance about zero is then 



a 2 (t) = 



1 



2t T 



t' 2 dt' = 



The error bar for the creation time becomes 



At = 



16. Ill 



ns = 23.9 ns. 



(2) 



(3) 



This is considerably larger than the final statistical error 
bar of 6.7 ns given in [l|. As seen from Figs. 4 and 11 of 
the reference, the explanation appears to be: The true 
proton distribution is not uniform and drops off at the 
beginning and end of the time interval. For our simple 
exercise we like the proton distribution to be flat, but let 
us use a new i max so that At becomes At = 6.7 ns. Solv- 
ing Q for i max gives i max = 1,473 ns, which is small on 
the scale of the figures, where 2 i max should cover about 
10,000 ns. Using t max = 5, 000 ns gives still 



At 



22.7 ns. 



(4) 



To explain the discrepancy with the error bar of the refer- 
ence to students, one has to take into account the struc- 
ture and statistical significance of the data and/or ana- 
lyze individual events at the boundaries. This requires 
access to the complete data sets underlying Fig. 11. 



Exercise 2: Let us give up the assumption that each 
proton creates neutrinos with the same probability. In- 
stead, let us assume that each proton diminishes by a 
factor 



< / = 1-p < 1 



(5) 



just a little bit the ability of the subsequent protons to 
create neutrinos. Let us assume that this is an uncorre- 
cted statistical process with the same / for each proton. 
What is the value of p, so that St of Eq. (JTJ is obtained? 

Under these assumptions the expectation value of t 
becomes 



(t) 



2t 



max 



te at dt 



2t 



max 



da 



dt 



d sinh(ai max ) 



da 



(6) 
(7) 
(8) 



For uniform creation probabilities a — and (t) = as 
it should. Now, we want to produce the effect ((TJ, which 
is small compared to i max . Therefore, the leading Taylor 
expansion of sinh(a£ max )/a is sufficient 



(*> = 



— e 
3 



da 



2 j-3 



a z t 



3! 



max q2 



(9) 



.(10) 



Solving the quadratic equation with £ max = 1, 473 ns and 
(t) = 60 ns gives two solutions for a of which the physical 
(discuss why) one is a = 9.67 x 10~ 5 ns -1 . 

As we assume our protons uniformly distributed and 
there are 6.2 x 10 15 protons in one extraction, the time 
interval associated with one proton is 



t P 



in 



2.37 x 10" 13 ns 



6.2 x 10 15 
and the factor / of Eq. ((SJ becomes 

/ = exp(-aip) = 1 

with 



P 



at P = 1.38 x 10~ 16 . 



(11) 



(12) 



(13) 



2 



Discussion: Is it obvious that < p <C 1.4 x 10 -16 
holds in the CNGS experiment? As described in detail 
by the CNGS collaboration [lj], the neutrino extraction 
succeeds in a cascade of events and the present authors 
understand nothing about the experimental details and 
whether any parts of the process could be affected by ini- 
tial protons. The purpose of our simple exercises can only 
be to provide a starting point for physics discussions with 
university undergraduates or students in gifted classes at 
high schools. Our idea is to address some of the issues of 
the experiment while staying within the calculus skills of 
our target group. Continuing our calculation within these 
limitations, the probability of the last proton to produce 
a neutrino is f n = exp(-a f max ) = 86.7% (n = 6.2 x 10 15 ) 
of the corresponding probability of the first proton. This 
number is visually too small when compared with the 
neutrino arrival time distribution of Fig. 11 of [l|. 

Students may want to vary the parameters used. For 
instance, if one takes i max = 5,000 ns as suggested by 
the quoted figures, assumes At = 22.5 ns and takes out 
two standard deviations (j4]) from the signal ([T]), one has 
only to explain St = 15 ns. Crunching the numbers again 
one ends up with an efficiency of /" = 99.1% for the last 
proton when compared to the first, which is no longer in 



visual discrepancy with Fig. 11. A possible modification 
of our approach is to let p depend on the proton num- 
ber i = 1, . . . , n, p(l) > p {p as calculated before) and 
p(i) — > with increasing i, a model suitable for computer 
simulations. 

Finally, one could ask the students to design an exper- 
iment that would efficiently rule out a bias in the depar- 
ture times. A solution would be to install a neutrino de- 
tector in the mountains close to CERN, because 8t does 
then not depend on the distance of the detector. Be- 
sides, critical students will presumably invent all kind of 
potential error sources, which can be ruled out by careful 
analysis. We hope that our exercises will serve educators 
and their students well to have a fun class or two about 
a hot topic. 
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